* 


Historic,  archived  document 

Do  not  assume  content  reflects  current 
scientific  knowledge,  policies,  or  practices. 


THE  INFLUENCE 
OF  A  FIRE -INDUCED 
CONVECTION  COLUMN 

ON  RADIOLOGICAL 
FALLOUT  PATTERNS 

A.  BROIDO 

A.  W.  McMASTERS 


U.S.D.A.  MAI 

AUG  2  1  ^OGo 
CATALOGING  PREP 


CALIFORNIA 
FOREST  AND  RANGE 
EXPERIMENT  STATION 

KEITH  ARNOLD,  DIRECTOR 
BERKELEY     -      MARCH  1959 


FOREST  SERVICE  -  U.  S.  DEPARTMENT  OF  AGRICULTURE 


ACKNOWLEDGMENTS 


This  study,  sponsored  "by  the  Federal  Civil  Defense  Administration 
(now  Office  of  Civil  and  Defense  Mobilization) ,  was  conducted  by  the 
Division  of  Forest  Fire  Research,  California  Forest  and  Range  Experiment 
Station,  under  subcontract  to  the  Institute  of  Engineering  Research, 
University  of  California,  Berkeley. 

The  authors  wish  to  acknowledge  the  assistance  of  Messrs.  William 
Go  O'Regan,  Robert  R.  Read,  and  Pinhas  Zusman  in  the  statistical  analyses 
and  computations  covered  so  briefly  in  this  report.    Other  participants 
in  various  phases  of  the  study  included  Wee  Yuey  Pong  and  Arthur  W. 
Jensen. 


THE  INFLUENCE  OF  A  FIRE- INDUCED  CONVECTION  COLUMN 
ON  RADIOLOGICAL  FALLOUT 


By  A.  Broido 
and  A.  W.  McMasters 

Division  of  Forest  Fire  Research 


The  material  used  in  this  report  was 
developed  under  a  research  contract 
sponsored  by  the  Office  of  Civil  and 
Defense  Mobilization. 


CALIFORNIA  FOREST  AND  RANGE  EXPERIMENT  STATION 
FOREST  SERVICE,  U.S.  DEPARTMENT  OF  AGRICULTURE 


The  Experiment  Station  is  maintained  at  Berkeley,  California 
in  cooperation  with  the  University  of  California. 


FOREWORD 


This  study  serves  as  another  illustration  of  how  the  diverse  re- 
search requirements  of  two  governmental  agencies  can  be  fulfilled  in  a 
single  research  program.    The  Forest  Service,  U.  S.  Department  of  Agricul- 
ture, as  a  result  of  its  responsibilities  in  the  control  of  forest  fires, 
has  considerable  interest  in  studying  the  wind  flows  in  and  around  fire- 
induced  convection  columns.    Such  flows  markedly  influence  subsequent 
fire  behavior.    One  way  to  study  the  wind  flows  aloft  is  through  the  in- 
troduction of  some  tracer  material  such  as  radiological  fallout.  Simul- 
taneously, the  Office  of  Civil  and  Defense  Mobilization  (OCDM)  has  an 
obvious  interest  in  the  effect  of  mass  fires  on  the  radiological  fallout 
distribution  pattern  to  be  expected  in  case  of  nuclear  attack.  Answers 
to  the  two  questions,  one  concerning  fallout  patterns  the  other  concern- 
ing wind  flows,  can  come  from  a  single  experiment.     Indeed,  from  such 
experiments  answers  to  either  question  would  be  meaningless  unless  the 
other  were  considered  also.    Thus,  continuation  of  the  program  begun 
with  this  study  under  OCDM  sponsorship  would  help  answer  some  of  the 
most  pressing  problems  in  radiological  defense  and  at  the  same  time  be 
of  application  to  those  in  the  Forest  Service  whose  responsibility  is 
prediction  of  fire  behavior  on  the  basis  of  wind  structure. 


ABSTRACT 


Since  no  nuclear  devices  have  been  detonated  "by  the  United  States 
under  conditions  leading  to  both  mass  fires  and  radiological  fallout,  a 
theoretical  and  small-scale  experimental  study  was  undertaken  to  see  if 
fire-induced  convection  columns  could  significantly  affect  fallout  pat- 
terns.   Experiments  were  conducted  in  a  6-  by  6-foot  low-velocity  wind 
tunnel  using  full-scale  fallout  simulants.    Values  of  characteristic 
parameters  of  fallout  patterns  were  compared  for  experiments  conducted 
with  and  without  fire,  all  other  factors  being  held  constant.    As  pre- 
dicted by  theory,  a  fire  upwind  of  the  region  of  deposition  introduced 
a  further  downwind  movement  and  lateral  dispersion  of  the  fallout.  In 
a  few  experiments  in  which  the  fire  was  located  at  the  edge  of  the  no- 
fire  pattern,  air  flow  into  the  convection  column  shifted  the  pattern 
peak  upwind.    These  results  suggest  quite  strongly  that  effects  which 
can  markedly  alter  civil  defense  planning  will  be  found  for  large-scale 
fires.    Since  this  effect  can  be  significantly  influenced  by  the  air  flow 
into  and  around  the  column,  open  air  experiments  on  a  larger  scale  should 
be  performed  so  that  such  flow  will  not  be  affected  by  tunnel  walls  and 
ceiling. 
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THE  INFLUENCE  OF  A  FIRE- INDUCED  CONVECTION  COLUMN 
ON  RADIOLOGICAL  FALLOUT  PATTERNS 


INTRODUCTION 

IMPORTANCE  OF  THE  STUDY 

The  detonation  of  a  high-yield  nuclear  weapon  at  or  near  the  sur- 
face of  the  earth  results  in  damage  which  has  been  summarized  in  terms  of 
four  principal  categories  of  effects:     (l)  blast,  (2)  thermal,  (3)  prompt 
radiation,  and  (4)  fallout.    The  first  two  may  act  jointly  in  producing  a 
mass  fire  which  can  extend  the  range  of  damage  well  beyond  the  extent  of 
the  primary  prompt  effects.    The  fourth  results  from  the  deposition  of 
large  quantities  of  highly  radioactive  particles  produced  in  the  detona- 
tion.   Distributed  by  air  currents,  these  particles  return  to  earth  as  a 
result  of  the  pull  of  gravity.    Radiation  from  this  fallout  extends  the 
range  of  damage  of  such  a  detonation  hundreds  of  miles  beyond  the  range 
of  the  other  three  effects  and  extends  the  time  of  significance  to  many 
days,  or  even  months,  after  the  time  of  detonation. 

Previous  studies  of  the  transport  of  small,  light  materials  such 
as  pollen  and  mold  spores  have  indicated  that  such  small  particles  are 
greatly  influenced  by  warm  air  currents  such  as  those  produced  over  an 
urban  complex.—/    Theoretical  analyses  of  convection  columns  above  large- 
scale  fires  indicate  that  updraft  velocities  resulting  from  even  moderate 
rates  of  heat  output  exceed  the  terminal  downward  velocities  of  most  radio- 
logical fallout  particles.     It  would  appear,  then,  that  convection  columns 
induced  by  fires  set  by  the  detonation  or  deliberately  set  after  the  deto- 
nation could  have  a  marked  effect  on  fallout  patterns.    Under  conditions 
of  sizeable  civil  defense  activity,  such  columns  could  greatly  alter  the 
necessary  civil  defense  planning  and  action. 

No  nuclear  detonations  by  the  United  States  have  permitted  obser- 
vation of  the  interrelationships  between  (l)  the  mass  fires  which  can 
result  from  such  detonations  and  (2)  the  fallout  of  radioactive  material 
following  a  surface  or  near-surface  burst.    Only  two  detonations  are  known 
to  have  occurred  under  conditions  where  mass  fires  could  form.    Both  these 
detonations,  at  Hiroshima  and  at  Nagasaki,  were  air  bursts  in  which  no 
fallout  would  be  expected.    All  test  devices  detonated  by  the  United 
States  have  been  set  off  at  sites  where  the  principal  surface  area  was 
sand  or  water,  neither  of  which  will  support  combustion.  Consequently, 
the  interrelationship  between  mass  fires  and  fallout  must  be  investi- 
gated theoretically  and  through  the  use  of  smaller-scale  experiments. 


1/    Heise,  H.  A.,  Meeting  of  the  Southern  Flying  Physicians 
Association,  Washington,  D.  C,  November  13,  1956.     (Reported  by  Asso- 
ciated Press.) 


PAST  WORK 


Harwit  (1956)  reported  the  results  of  a  study  on  the  deposition 
of  airborne  particles  on  a  large  heated  surface.    This  study  considered 
cases  of  low  heat  output  (of  the  order  of  micro  calories  per  sq.  cm.  per 
sec;  for  example,  the  natural  heating  of  an  urban  area  to  a  few  degrees 
above  the  ambient  temperature  of  the  surroundings) .    From  theoretical 
computations ,  Harwit  estimated  that  heat  outputs  greater  by  at  least  a 
factor  of  a  hundred,  (producing  upward  velocities  of  the  order  of  100  cm. 
per  sec.  or  more),  would  be  needed  to  produce  an  appreciable  effect  on 
the  particle  deposition  pattern. 

Theoretical  treatment  of  the  convection  above  heat  sources  (Scesa 
&  Sauer,  1954;  Scesa,  1957)  indicate,  too,  that  when  the  heat  output  of 
the  source  is  low  with  respect  to  the  ambient  winds,  the  effect  of  buoy- 
ancy may  be  neglected  (no  convection  column  will  form) >    On  the  other 
hand,  heat  outputs  which  may  be  expected  in  fire  situations  similar  to 
those  in  forest  fires  (of  the  order  of  kilo-calories  per  sq.  cm.  per 
sec.)  are  sufficient  to  permit  the  formation  of  a  convection  column  in 
spite  of  the  presence  of  moderately  strong  natural  winds.     In  such  con- 
vection columns  the  buoyancy  effect  produces  updraft  velocities  of  the 
order  of  thousands  of  cm.  per  sec.    Since  the  terminal  velocities  of  most 
fallout  particles  are  expected  to  be  lower  than  these  values  (Kellogg,  et 
al.,  1956),  such  columns  should  have  a  marked  influence  on  the  deposition 
pattern  of  such  fallout  particles „    For  fires  started  by  the  nuclear  deto- 
nation itself,  such  influence  may  include: 

1.  Prevention  of  deposition  in  an  area  which  might  otherwise 
be  expected  to  be  highly  radioactive. 

2.  Concentration  of  fallout  in  regions  which  might  otherwise 
be  relatively  less  radioactive. 

3.  Under  certain  circumstances,  a  gross  change  in  the  entire 
fallout  pattern. 

In  addition,  certain  situations  might  lend  themselves  to  the  use  of  delib- 
erately set  fires  as  contamination  prevention  or  decontamination  tools. 

In  September  of  1956,  a  series  of  preliminary  experiments  were  per- 
formed to  investigate  the  possibility  of  burning  out  the  fuel  to  decon- 
taminate a  combustible  area  which  had  been  exposed  to  radioactive  fallout 
(Broido,  1957)°    Although  the  experiments  were  limited,  the  broad  con- 
clusions reached  included  the  following:    In  areas  with  fuel  content  below 
1  ton  per  acre,  only  a  small  degree  of  decontamination  (less  than  20%) 
could  be  expected.    Fast-forming,  rapidly  burning  fires  in  contaminated 
areas  as  small  as  1/10  of  an  acre  with  fuel  loading  of  20  tons  per  acre 
can  result  in  the  removal  of  perhaps  a  third  of  the  fallout,  and  the  re- 
moved material  is  disbursed  so  that  there  is  no  significant  concentration 
in  any  other  region.    The  experiments  indicated  that  under  the  much  more 
severe  conditions  found  in  a  large  area  mass  fire,  more  striking  effects 
might  well  occur. 
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Shipley  (1958)  published  the  results  of  a  brief  study  on  the  dis- 
persal of  fallout  by  means  of  a  heat  source .    His  theoretical  analysis 
included  a  procedure  which  admittedly  was  "arbitrary  and  had  no  theoretical 
or  experimental  basis  for  being  correct."    His  experimental  technique  was 
to  measure  the  number  of  particles  of  chalk  dust  falling  on  slides  set 
under  and  near  an  electrically  heated  wire  grid  mounted  in  a  small  office 
room.    In  these  experiments  the  chalk  dust  was  observed  to  be  lifted  as  it 
passed  over  the  heat  source  area.    However,  a  count  of  the  fallout  deposited 
on  and  around  the  area  revealed  only  a  small  decrease  in  the  amount  of 
material  deposited.    The  author  points  out  that  these  unfavorable  results 
may  have  been  the  fault  of  confining  the  experiment  to  a  small  room.  He 
states:     "Had  the  experiments  been  conducted  in  a  much  larger  area  it  may 
be  that  fallout  dissipated  from  the  heat  sources  could  have  been  carried 
further  away  from  the  protected  area  and  there  deposited." 

SCOPE  OF  THE  PRESENT  STUDY 


A  preliminary  analytical  evaluation  of  the  problem  indicated  that 
no  true  small-scale  modelling  was  feasible.    However,  the  availability  of 
a  6-  by  6- foot  low  velocity  wind  tunnel  which  permitted  the  introduction 
of  fires  giving  appropriate  updraft  velocities  made  it  appear  that  full- 
scale  particles  having  distribution  as  actually  observed  in  fallout  might 
give  results  of  considerable  interest. 


THEORETICAL 


EFFECT  OF  AN  IDEALIZED  CONVECTION  COLUMN 
ON  THE  FALL  OF  A  SINGLE  PARTICLE 


Assume  a  particle  at  point    A    in  figure  1,  a  distance    H  above 
ground  in  an  idealized  homogeneous  atmosphere  of  constant  density1,  falling 
with  a  terminal  velocity    V     in  a  wind  of  constant  velocity    V  .    The  vec 
tor  diagram  is  simply:  ^ 


and  the  particle  would  fall  along  the  path    AR    in  figure  1,  landing  at 
point    R  ,  a  d: 
H  /V    ,  and  so 


point    R  ,  a  distance    Lq    downwind  of  point  A.    The  time  of  fall  would  be 


V  H 

L    =  V  t  =  -ip  (1) 

o       w  V 
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Figure  l.--Fall  of  a  single  particle  through  an  idealized 
convection  column. 
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Now,  assume  an  idealized  convection  column  is  introduced  into  the 
system.    This  convection  column  is  produced  by  a  source  generating  air  of 
the  same  physical  properties  as  the  surroundings  (temperature,  pressure, 
etc.)  and  pumping  this  air  upwards  at  a  velocity    Vc  .    This  column  of  air 
rises  at  an  angle  dependent  upon    Vc    and    Vw    and  is  then  fed  into  a  sink. 
A  slip  surface  is  assumed  to  exist  at  the  column  boundary,  so  that  the  pre- 
sence of  the  column  introduces  no  circulation  in  the  existing  atmosphere,, 
Thus,  the  trajectory  of  the  particle  would  not  change  until  it  reached  the 
leading  edge  of  the  column  (at  a  height    H,  ) .    At  this  point  the  vertical 
velocity  of  the  particle  with  respect  to  the  ground  becomes    Vp  -  Vc  . 
If    Vp  >  Vc  ,  the  particle  would  continue  to  fall,  but  more  slowly.  If 
Vp  =  vc  ,  the  particle  would  move  through  the  column  in  a  horizontal  path 
with  velocity    Vw  .     If    Vp  <  V0  ,  the  particle  would  rise  as  it  moved 
through  the  column  in  a  path  determined  by  the  vector  diagram: 


V 
w 


V  -V 
c  p 


and  shown  as  the  dashed  line  within  the  column  in  figure  1.    When  the  par 
ticle  reached  the  downwind  edge  of  the  column,  it  would  again  fall  at  the 
same  angle  as  before  (following  the  dashed  line  in  figure  l)  and  land  at 
point    S  .    The  downwind  distance  of  fall,    L  ,  can  be  obtained  by  com- 
puting   hi  j  L2  t  and  L-^    and  summing  them. 

As  before,    L-^    and    L-^    are  given  simply  as 

V  (H  -H  )  V  H 

,  W        O         1  n         t  W  d 

Ll  =  V   and    L3  =  ~V~~ 

P  P 

In  a  similar  manner,  two  expressions  may  be  found  for  : 

L2  =  —   and    L2  =  Lo  +   V  

c    p  c 


The  displacement  of  the  particles  caused  by  the  convection  column, 


AL  =  L 


may  be  obtained  from  the  above  equations  as 


AL  =  L1  +  L2  +  L3 


fV 

L    =  L  ItT 

o       c  V 

P 


(2) 
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Thus,  the  displacement  of  the  particle  by  such  an  idealized  con- 
vection column  is  seen  to  be  completely  independent  of  the  wind  velocity, 
being  simply  a  function  of  the  column  width  and  updraft  velocity  and  the 
fall  velocity  of  the  particle. 

If  the  physical  characteristics  of  the  column  are  allowed  to  change 
(e.g.,  heated  air),  the  fall  velocity,  with  respect  to  the  column  gases, 
of  the  particle  within  the  column  is  not  the  same  as  that  outside  the 
column.     In  this  case,  the  displacement  of  the  particle  becomes: 


AL  =  L 


v2  +  V  (V  -V' ) 
c        c    p  p' 

V  V' 
P  P 


(3) 


where  V"  is  the  terminal  velocity  of  the  particle  in  the  medium  of  the 
column.  ^ 


THE  EFFECT  OF  AN  IDEALIZED  CONVECTION  COLUMN  ON 
AN  EXTENDED  SOURCE  OF  FALLOUT  PARTICLES 


Line  Source  of  Particles 


Consider  the  fallout  particles  in  the  previous  example  to  be  one 
of  a  large  number  of  identical  particles  lined  up  in  the  same  horizontal 
plane  (fig.  2).     In  falling  from  height    HQ  to  Hi    these  particles  would 
follow  parallel  paths  uninfluenced  by  the  existence  of  the  convection 
column  shown  in  figure  2.    At  height    Hi    the  particle  at    Ax  enters  the 
convection  column  and,  following  the  path  indicated  by  the  broken  line, 
ultimately  lands  at  point    S     (rather  than  at  point    R    where  it  would 
land  had  no  convection  column  been  present).    The  following  particles 
enter  the  column  at  progressively  later  times  and  at  lower  elevations. 
Upon  entering  the  column  they  follow  a  path  parallel  to  that  taken  by 
the  first  particle  and  their  points  of  deposition  are  displaced  accord- 
ingly.   Thus,  the  particle  at    P/  when  the  first  particle  entered  the 
column  enters  upon  reaching  point    M    and  is  displaced  to  land  at 
point    R  .    A  particle  starting  its  fall  at  point    C    would  enter  the 
column  at  its  base  (point  P)  and  like  the  previous  particles  be  dis- 
placed, so  as  to  land  at  point    Q  .     The  particle  at    D  ,  however,  and 
all  subsequent  particles,  would  follow  a  straight- line  path  directly  to 
the  ground- landing  in  exactly  the  same  position  whether  or  not  the  column 
were  present. 

Thus,  without  a  column  the  line  of  particles,  EA  ,  would  land  in 
a  line  of  equal  length,  NR  ,  on  the  ground.  The  presence  of  the  column 
does  not  in  any  way  change  the  length  of  this  line,  but  it  does  break  it 
into  two  segments,     NO    and    QS  . 

The  effect  of  the  convection  column  on  the  accumulation  of  fall- 
out particles  on  the  ground  in  such  a  case,  then,  is  dependent  upon  the 
point  under  consideration.    Between  points    N    and    0    no  change  occurs, 
that  is,  the  same  particles  land  at  exactly  the  same  points.     From  P 


-6- 


Figure  2. --Fall  of  an  extended  source  of  identical  particles. 
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to    Q    fallout  particles  which  would  arrive  without  a  convection  column 
are  totally  displaced  and  no  particles  land  in  this  region .    Prom    Q  to 
R  ,  although  a  displacement  occurs,  the  particle  which  would  normally 
land  in  this  region  is  replaced  by  another  identical  particle  and,  thus, 
no  change  is  observed..    From    R    to    S  ,  particles  now  land  where,  with- 
out a  convection  column,  no  particles  would  have  been  found. 

The  previous  paragraph  illustrates  the  effect  of  such  an  idealized 
convection  column  on  a  line  of  particles  which  are  identical  in  size  and 
fall  velocity.    Consider,  now^  another  such  line  of  smaller  particles 
located  at  such  a  position  that,  with  their  particular  fall  trajectory, 
they,  too,  would  land  in  the  region    NR    if  no  column  were  present.  The 
effect  of  the  convection  column  would  be  qualitatively  the  same  as  that 
given  in  the  previous  paragraph,  but  these  smaller  particles  would  be  dis- 
placed a  greater  distance  by  the  column.    The  degree  of  overlap  of  the 
downwind  segments  would  be  a  function  of  the  terminal  velocity  of  the  con- 
vection column.    If,  then,  a  large  number  of  different  sized  particles 
were  to  be  displaced  by  such  an  idealized  column,  the  over-all  result  in 
the  downwind  region  would  be  an  increase  in  the  size  of  the  region  on 
which  fallout  would  land,  with  a  corresponding  decrease  in  concentration 
at  any  one  point. 

Cloud  of  Particles 

The  extension  of  the  argument  to  a  3 -dimensional  cloud  of  particles 
is  straightforward.    Consider,  for  example,  the  area    EAA'e'     in  figure  2 
to  represent  a  cross-sectional  slice  through  an  idealized  cloud  of  identi- 
cal particles.    Without  a  convection  column  all  of  the  particles  in  the 
line    M'  would  land  at  point    R  .    Those  in  the  line    EE1  would  land  at 
point    N  .    Under  the  influence  of  the  idealized  convection  column,  all 
of  the  particles  in  the  line    AA!  now  land  at  point    S  .    However,  the 
particles  in  the  line    BB'  now  land  at  point    R    in  place  of  those  which 
have  been  displaced.    Thus,  assuming  a  uniform  concentration  through  the 
cloud,  if  the  lines    BB1  and    AA!  are  equal  in  length,  no  change  will 
occur  in  the  concentration  of  particles  landing  at  point    R  .  Similarly, 
for  this  example,  the  effects  at  any  other  point  on  the  ground  follow 
directly  from  what  has  already  been  said  about  the  line  source. 

This  example,  however,  has  considerably  idealized  the  true  situation 
by  assuming  the  cloud  to  have  constant  horizontal  upper  and  lower  boundaries 
and  side  boundaries  with  slopes  identical  to  the  fall  trajectory  of  the  par- 
ticles.   Consideration  of  the  real  dimensions  of  a  fallout  cloud  as  well  as 
expected  wind  velocities  and  fall  velocities  of  the  particles  would  indi- 
cate that  a  more  reasonable  picture  is  that  given  in  figure  3  °    Tne  wind 
velocities  are,  in  general,  expected  to  be  considerably  greater  than  the 
fall  velocities  of  the  particle,,  and  thus  the  trajectory  would  be  more 
nearly  horizontal.    In  such  cases,  the  line  through  the  cloud  of  particles 
which  would  land  at  a  given  point  on  the  ground  would  terminate  at  the  up- 
wind edge  rather  than  the  upper  edge  of  the  cloud a    A  downwind  shift  of 
this  cloud,  then,  would  result  in  a  shortening  of  such  a  line  and  a  de- 
crease in  concentration  of  particles  at  a  point.    Thus,  the  introduction 
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of  the  convection  column  shown  in  figure  3  would  result  in  the  particles 
in  the  line    BB'  landing  at  the  point  which  otherwise  would  be  occupied 
by  those  in  the  line    AA1   .    Since    BB'  is  just  half  the  length  of    AA'  , 
the  concentration  of  particles  of  this  size  landing  at  point    R    would  be 
cut  in  half .2/ 

THE  EFFECT  OF  A  REAL  CONVECTION  COLUMN 

The  factors  which  have  been  considered  up  to  this  point  are  those 
which  influence  the  fall  of  the  particles  within  the  convection  column  it- 
self,,   In  the  idealized  example  considered,  no  other  effects  could  occur 
since  the  gases  of  the  convection  column  were  generated  and  dispelled 
without  being  circulated  into  the  surrounding  atmosphere.    In  a  real  con- 
vection column,  such  circulation  must  be  taken  into  account. 

Although  some  gases  are  produced  in  a  combustion  reaction,  much  of 
the  gas  rising  in  a  fire-induced  convection  column  is  air  entering  the 
column  at  various  elevations  above  grade.    The  rising  gases  do  not  termi- 
nate at  some  sink,  but  recirculate  and  influence  the  general  flow  of  air 
in  a  large  area  around  the  column  itself,    Corcos  (1958)  has  shown  that 
particles  of  the  type  considered  here  follow  quite  closely  the  circulation 
patterns  of  the  atmosphere.    Thus,  before  the  total  effect  of  a  convection 
column  on  a  fallout  pattern  can  be  ascertained,  the  flow  of  air  around  the 
convection  column  must  be  understood. 

Unfortunately,  practically  no  information  is  available  about  such 
flow.    Strong  indrafts  have  been  measured  at  the  base  of  a  fire-induced 
column,  but  little  is  known  about  such  questions  as:    What  fraction  of  the 
air  entering  the  column  enters  at  the  base  and  what  fraction  at  higher  ele- 
vations?   What  fraction  enters  downward  of  the  fire  and  what  fraction  up- 
wind or  crosswind?    When  the  rising  air  "spills  over"  at  the  top  of  the 
column,  how  far  does  it  travel  at  these  upper  elevations  before  recycling 
at  elevations  near  the  ground?    Until  more  is  learned  about  answers  to 
questions  of  this  type,  little  more  can  be  ascertained  theoretically  about 
the  influence  of  a  real  convection  column  on  a  fallout  pattern.    In  fact, 
it  is  to  be  hoped  that  experimental  studies  of  such  influences  may  help 
provide  additional  insight  into  these  air  flows, 

EXPERIMENTAL 
FALLOUT  SIMULANT 

When  a  nuclear  weapon  is  detonated  on  or  near  the  surface,  great 
masses  of  pulverized  debris  are  sucked  upward  from  the  surface  and  are 
mixed  with  radioactive  fission  products  released  by  the  explosion.  The 
composition  of  this  fallout  depends,  of  course,  upon  the  composition  of 
the  surface  over  which  the  detonation  took  place.    The  fallout  particles 
vary  considerably  in  shape  and  in  size.    Consequently,  the  rate  of  de- 
scent of  fallout  ranges  from  the  "throw-out,"  the  rapid  deposition  of 

2/    A  somewhat  different  treatment  is  given  by  Corcos  (1958),  who 
considers  the  concentration  gradient  of  particles  in  a  non-uniform  cloud. 
Application  of  this  treatment  to  the  discontinuous  cloud  considered  here 
would  give  identical  results  if  the  concentration  gradient  were  taken  as 
that  over  a  line  the  length    AA'  extended  through  the  points    BB1  and  on 
beyond  the  cloud. 
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large-sized  material  in  the  vicinity  of  the  crater,  to  the  slow  rate  of 
descent  of  the  smallest  particles  which  circle  the  earth  many  times. 
This  study  dealt  with  the  intermediate  range  of  particle  sizes,  namely, 
those  that  would  he  expected  "beyond  the  region  of  total  primary  destruc- 
tion, hut  whose  concentration  would  be  sufficiently  high  to  create  a 
severe  radiological  hazard.    Specifically,  attention  was  focused  on  par- 
ticles in  the  100-  to  300-micron  range. 

Choice  of  Simulant 

Three  principal  factors  influenced  the  selection  of  a  simulant  to 
this  fallout  material: 

1.  Theoretical  computation  (see,  for  example,  Corcos, 
(1958),  indicated  that  a  true  scaling  experiment 
was  not  feasible. 

2.  The  use  of  very  much  lighter  particles  would  intro- 
duce the  possibility  of  other  complicating  effects 
and  make  extrapolation  more  questionable. 

3.  The  U.  S.  Naval  Radiological  Defense  Laboratory 
(U.S.N.R.D.L. )  was  using,  and  had  available,  a 
simulant  of  appropriate  composition,  density, 
and  shape  and  size  distribution.    This  material, 
also  used  in  the  earlier  decontamination  experi- 
ments (Broido,  1957),  consisted  of  soil  taken 
from  the  test  site  at  Camp  Stoneman,  California 
and  dried,  pulverized,  and  sieved  until  it  passed 
a  No.  30  mesh  screen. 

To  simplify  interpretation  of  the  results,  this  material  was  further 
subdivided  by  sieving  through  a  series  of  Tyler  Standard  Screens.    In  this 
procedure,  very  little  material  was  found  of  size  large  enough  to  be  held 
by  a  250-|_l  screen,  and  the  sieving  procedure  became  too  slow  to  be  practi- 
cal for  sieve  sizes  smaller  than  74-M..    Since  this  smaller  material  settled 
too  slowly  to  be  used  within  the  confines  of  the  wind  tunnel  available,  no 
further  breakdown  was  attempted. 

The  fallout  material  used,  then,  consisted  of  five  fractions  of 
material  not  held  by  a  250-(j.  screen  and  held,  respectively,  by  210-u,  177-H> 
149-i-L,  105-i-i,  and  T4--M-  screens. 

Method  of  Measuring 

An  appropriate  method  for  determining,  or  even  defining,  the  "size"  of 
an  irregularly- shaped  particle  is  subject  to  considerable  arbritrariness . 
Thus,  a  sieve  with  square  holes  of  side         in  length  will  pass  spherical 
particles  of  diameter  up  to  length  On  the  other  hand,  a  particle  in 

the  form  of  a  square  plate  can,  in  principle,  pass  through  such  a  hole  if 
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the  length  of  the  edge  of  the  plate  is  less  than  V2"  JL .    Measured  under  a 
microscope,  then,  at  least  some  dimensions  of  an  irregularly- shaped  parti- 
cle passing  through  a  250-u  screen  would  he  larger  than  250-u. 

Even  if  all  of  the  dimensions  of  an  irregularly-shaped  particle  were 
accurately  known,  the  definition  of  the  particle  "diameter"  would  be  sub- 
ject to  considerable  ambiguity  (Cadle,  1955) »    Statistical  diameters,  based 
upon  measuring  the  dimensions  of  a  large  number  of  particles,  are  usually 
used.    If  the  length,  breadth,  and  thickness  of  the  particle  can  be  esti- 
mated, the  arithmetic  average  of  the  three  diameters  may  be  used.  Alter- 
natively, the  diameter  of  a  sphere  (or  the  length  of  the  side  of  a  cube) 
of  equal  volume  may  be  defined  as  the  particle  diameter.    A  useful  defini- 
tion of  diameter  of  particles  measured  under  a  microscope  was  proposed  by 
Martin  (1924) .     It  is  the  distance  between  opposite  sides  of  the  particle 
measured  crosswise  of  the  particles  and  on  a  line  bisecting  the  projected 
area.    Since  irregularly- shaped  particles  will  generally  land  "flat," 
that  is,  with  their  smallest  dimensions  perpendicular  to  the  surface  on 
which  they  land,  this  measurement  will  not  give  a  good  indication  of  the 
smallest  dimension.    However,  if  a  sufficiently  large  number  of  particles 
are  measured,  a  good  average  value  of  the  other  two  dimensions  of  the 
particle  will  result. 

In  the  course  of  these  experiments  selected  samples  were  collected 
on  gummed  transparent  paper.    After  a  preliminary  statistical  analysis 
established  that  these  particles  were  randomly  distributed  throughout  the 
paper,  mean  Martin's  diameters  were  obtained  by  dividing  the  paper  into  a 
large  number  of  small  squares  and  selecting  the  squares  for  measurement 
from  a  table  of  random  numbers.    These  diameters  were  then  compared  with 
the  gross  values  obtained  using  the  sieving  technique. 

Since,  for  purposes  of  these  experiments,  the  particle  property  of 
prime  interest  is  the  fall  velocity  of  the  particle,  an  experiment  was  per- 
formed to  establish  an  "effective  diameter."    A  vertical  wind  tunnel,  to 
be  described  in  a  subsequent  report,  was  designed  and  built  (fig.  4) • 
Fallout  particles  were  then  mixed  with  spherical  glass  beads  (fig.  5)  of 
appropriate  size.    The  mixture  was  fractionated  in  the  vertical  wind  tun- 
nel and  the  various  fractions  measured  under  the  microscope  (fig.  6). 
The  diameter  of  the  glass  beads  in  a  given  fraction  established  the  fall 
velocity  of  the  entire  fraction.    The  Martin's  diameter  of  the  fallout 
could  then  be  related  to  this  effective  diameter . 

Table  1  gives  some  representative  results  of  such  a  separation. 
For  these  experiments,  a  gross  sieve-size  sample  was  fed  into  the  hopper 
of  the  vertical  wind  tunnel.    Each  such  sample  was  then  fractionated  by 
progressively  lowering  the  voltage  to  the  wind  tunnel  exhaust  blower  in 
5  volt  increments  and  collecting  the  material  falling  out  of  the  tunnel. 
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Exhaust  fan 


Sample  hopper- 


Figure  4. --Vertical  wind  tunnel 
used  to  fractionate  irregularly 
shaped  particles  in  terms  of 
fall  velocity. 


Contraction 
section  


Sample  falling  //  / 


Diffusion 
section 
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Figure  5. --Mixture  of  glass  beads  added  to 
fallout  particles  prior  to  separation. 


Figure  6. --Typical  sample  of  fractionated 
material.     Constant  size  of  glass  beads 
establishes  fall  velocity  of  irregularly 
shaped  fallout  particles. 
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Method  of  Dispersing 


Mathematical  fallout  models  used  to  predict  the  distribution  of 
fallout  on  the  ground  after  a  nuclear  detonation  generally  begin  with  a 
fallout  cloud  stabilized  at  its  maximum  elevation . '  The  models  then  at- 
tempt to  predict  the  motion  of  the  particles  constituting  this  cloud  as 
a  result  of  the  pull  of  gravity  and  the  influence  of  the  wind.  Such 
models  eliminate  the  complicating  features  introduced  by  the  detonation 
itself  and  by  other  parameters  important  during  the  formative  stages  of 
the  cloud.    Similar ly,  these  experimental  studies  attempted  to  produce 
a  cloud  which  would  have  stabilized  and  which  would  be  moving  only  as  a 
result  of  the  influences  of  gravity  and  the  wind  at  the  time  the  parti- 
cles entered  a  convection  column.    The  ideal  experimental  cloud  would  be 
completely  uniform,  particularly  in  its  crosswind  dimension,  and  would 
be  wider  than  the  convection  column  used.    In  addition,  the  concentration 
of  particles  would  not  be  sufficient  to  influence  the  physical  character- 
istics of  the  air  through  which  they  fell,  but  would  be  sufficiently  great 
to  give  a  readily  measurable  pattern., 

The  simplest  method  of  injecting  the  fallout  particles  into  the  wind 
tunnel  airstream  would  be  to  drop  the  sample  from  a  funnel.    This,  however, 
would  produce  an  undesirably  narrow  pattern.    Dropping  from  a  hopper  orient- 
ed lengthwise  across  the  tunnel  was  ruled  out  because  of  the  inability  of 
getting  a  uniform  flow  from  such  a  hopper.    The  various  devices  used  by  the 
U.S.N.R.D.L.  to  spread  contaminant  in  their  decontamination  studies  were 
examined,  and  the  only  one  which  showed  promise  for  these  experiments,  a 
De  Vilbiss  Powder  Blower,  Model  175,  was  tried.    However,  this  device  proved 
to  be  too  small  to  produce  a  sufficiently  large  cloud  of  material. 

The  disperser  finally  selected  was  a  Binks  Model  No.  171  Flock  Gun 
(fig.  7) •    This  compressed-air-powered  spray  gun  was  designed  to  handle  par- 
ticulate matter  mixed  with  paints.    The  metal  can  that  came  with  the  gun 
was  removed  and  the  gun  was  adapted  to  operate  with  the  sample  contained  in 
a  one-quart  glass  jar.    A  series  of  experiments  were  performed  to  investi- 
gate the  effect  of  air  pressure,  level  of  fallout  in  the  jar,  gun  nozzle 
setting,  and  duration  of  firing  on  the  shape  and  consistency  of  the  pattern. 
Examination  of  the  data  indicated  that  the  most  satisfactory  results  were 
obtained  with  the  jar  filled  to  one-third  to  two-thirds  capacity,  with  the 
air  pressure  at  30  psig  for  a  10-second  run,  and  with  the  gun  nozzle  adjust- 
ed to  give  a  wide,  flat  horizontal  cloud. 

Although  the  carefully  controlled  conditions  led  to  a  fairly  repro- 
ducible fallout  pattern,  the  total  quantity  of  material  discharged  from  the 
gun  in  a  run  of  fixed  time  duration  varied  by  as  much  as  a  factor  of  5  for 
successive  runs.    Thus,  for  all  experiments,  the  general  shape  of  the  dis- 
tribution, rather  than  the  total  amount  of  material  leaving  the  gun,  was  of 
significance. 
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Since  the  flock  gun  was  to  "be  used  in  a  6-foot  wind  tunnel,  the  maxi- 
mum height  above  ground  at  which  it  could  be  mounted  was  about  5-1/2  feet. 
Fired  under  standard  conditions  from  this  height,  the  gun  gave  a  pattern 
that  was  roughly  2  to  3  feet  wide  and  about  10  feet  long.    The  lengthwise 
distribution  of  material  in  this  pattern  was  log  normal.    Thus,  the  charac- 
teristic parameters  of  the  distribution  patterns  which  could  be  measured 
with  fair  consistency  were  the  positions  of  the  logarithmic  mean,  the  vari- 
ance, and  the  peak  in  the  distribution  curve. 

A  desirable  characteristic  of  a  satisfactory  fallout  disperser  is 
that  it  not  influence  the  results  of  the  wind  tunnel  experiment  to  be  per- 
formed.   This  means  that  the  dispersive  effects  of  the  flock  gun  should 
dissipate  rapidly  and  the  cloud  of  fallout  should  be  falling  vertically, 
free  from  the  influence  of  the  gun,  well  before  it  hits  the  ground.    As  a 
check,  a  series  of  experiments  were  performed  in  still  air  with  the  various 
particle  sizes  and  with  the  gun  mounted  at  distances  from  9  inches  to  5-1/2 
feet  above  a  table  top.    Cross-pattern  strips  of  paper  3  inches  wide  were 
mounted  on  the  top  of  the  table,  and  samples  were  collected  from  these 
strips  and  weighed.    The  cumulative  weights  were  plotted  on  log  probability 
paper  and  the  position  of  the  pattern  peak  then  plotted  against  height  of 
disperser  above  table  top. 

The  gun  effect  dissipated  quite  rapidly  and  particles  were  essential- 
ly in  free-fall  after  they  had  fallen  only  a  foot  or  so  (fig.  8).  Even 
more  favorable  results  would  be  expected  under  conditions  of  horizontal  wind 
since  upon  leaving  the  gun  the  particle  velocity  would  need  decrease  only 
to  the  wind  velocity  through  the  tunnel  and  not  to  zero  velocity  as  it  must 
in  still  air. 

HEAT  SOURCE 

These  experiments  required  a  heat  source  that  would  produce  a  con- 
vection column  2  or  3  feet  across  with  updraft  velocities  of  the  order  of 
several  miles  per  hour.    Also,  the  flames  had  to  be  readily  lighted  and 
extinguished  to  minimize  delays  between  experiments.    Gas  flames  meet  this 
latter  requirement  ideally.    Since  cylinders  of  propane  gas  were  readily 
available  at  the  wind  tunnel  site,  this  gas  was  chosen  as  the  fuel  for  all 
fire  experiments.    The  heat  source  consisted  of  a  group  of  four  National 
Welding  Equipment  Company  ribbon  flame  burners.    Each  had  a  burning  area 
3 A  inch  by  12  inches,  and  the  burning  surface  was  located  9-1/2  inches 
above  the  floor.    The  original  plan  was  to  conduct  the  experiments  with 
two  rows  of  burners  positioned  with  their  axes  at  right  angles  to  the  wind 
tunnel.    However,  it  was  found  that  in  this  position  the  stronger  winds 
readily  extinguished  the  flames.    Consequently,  the  basic  configuration 
consisted  of  two  rows  of  burners  with  axes  parallel  to  the  axis  of  the  tun- 
nel and  positioned  as  close  together  as  the  burner  bases  would  allow,  about 
10  inches  between  burners  (fig.  9). 
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LOCATION   OF   PEAK  (feet) 


Figure  8. --Particle  trajectory  when  fired  from  fallout  disperser 
into  still  air. 
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Propane  gas  leaving  the  cylinder  was  passed  through  a  Flamo  pres- 
sure regulator  which  held  the  output  pressure  at  4  p.s.i.    The  gas  was 
then  fed  through  1/2-inch  inside  diameter  rubber  hose  to  the  burners.  A 
flowmeter  placed  in  the  gas  line  gave  the  flow  rate  of  gas  as  I.78  cubic 
feet  per  minute.    Assuming  complete  combustion,  this  corresponded  to  a 
heat  output  of  4,500  B.T.U.  per  minute„    Although  the  burners  had  facil- 
ities for  premixing  of  the  gas  with  compressed  air,  no  such  premixing 
was  used  since  more  satisfactory  characteristics  of  the  convection  column 
were  obtained  under  free-burning  conditions. 

Under  the  conditions  of  use  the  gas  burned  with  a  yellow  flame 
about  2  feet  long.    The  parallel  flames  from  the  two  sets  of  burners  seem- 
ed to  blend  into  one  large  flame  and  a  single  convection  column.  However, 
as  a  check  against  a  possible  valley  in  the  center  of  the  column,  some 
experiments  were  performed  with  the  front  burners  tilted  toward  each  other 
in  the  shape  of  a  "V"  pointing  in  the  direction  in  which  the  wind  was  blow- 
ing (fig.  10). 

Under  both  configurations  the  cross-sectional  area  of  the  flames  was 
approximately  2  square  feet.    For  this  area  the  heat  output  of  the  burners 
corresponded  to  the  burning  of  wood  fuels  at  the  rate  of  approximately 
6  tons  per  acre  per  minute. 

Preliminary  estimates  were  made  of  the  updraft  velocities  of  the 
convection  columns  in  the  wind  tunnel  with  a  wind  velocity  of  6  miles  per 
hour .    Measurements  made  with  a  Byram  vane  anemometer  held  with  its  axis 
in  a  vertical  position  at  a  point  8  feet  downstream  from  the  burners  and 
4-1/2  feet  above  the  tunnel  floor  indicated  an  updraft  velocity  of  about 
2.2  miles  per  hour.    This  measurement  is  confirmed  roughly  also  by  the 
angle  the  flame  made  with  the  horizontal  wind  stream  and  by  anemometer 
measurements  taken  at  this  angle.    Unfortunately  no  instrumentation  is 
available  to  permit  accurate  measurement  of  wind  velocities  in  hot  con- 
vection columns.    Further  work  is  being  done  to  develop  an  instrument  for 
measuring  these  velocities  more  accurately. 

WIND  TUNNEL  EXPERIMENTS 

Description  of  Tunnel 

The  wind  tunnel  experiments  were  conducted  in  a  low- velocity  Eiffel- 
type  tunnel  located  outdoors  at  the  San  Dimas  Experimental  Forest  in  southern 
California  (fig.  11).    The  tunnel,  described  in  detail  by  Fons  (1940)  is 
52  feet  long  and  is  composed  of  three  main  segments.    The  bellmouth  entrance 
is  5  feet  long,  converging  from  11-1/2  feet  square  to  6  feet  square.  The 
test  section  is  6  feet  square  and  32  feet  long  (an  8-foot  extension  of  the 
test  section  was  introduced  since  its  description  by  Fons)„    To  help  insure 
parallel  flow  of  air  in  the  test  section^  a  honeycomb  is  placed  at  the  test 
section  entrance.    The  exit  consists  of  a  converging  section,  a  honeycomb 
and  screen,  and  a  diverging  section.    A  four-bladed  propeller-type  fan, 
4  feet  in  diameter,  is  located  at  the  exit  of  the  diverging  section.  The 
fan  is  belt-driven  by  a  six-cylinder  automobile  engine,  with  a  tachometer 
mounted  on  the  propeller  shaft  to  determine  the  fan  speed. 
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Figure  10. --"V"  configuration  of  gas  burners 
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Before  these  experiments  the  wind  tunnel  had  last  "been  used  in  the 
spring  of  1945.    Despite  its  exposure  to  the  elements  and  minimal  mainten- 
ance ,  all  of  the  equipment,  including  the  automobile  engine  used  to  drive 
the  fan,  was  still  in  operating  condition.    A  limited  number  of  calibration 
checks  of  wind  speed  versus  propeller  speed  in  the  2-  to  6-mile  per  hour 
range  indicated  that  previous  calibrations  could  be  applied  without  correc- 
tion.   In  fact,  only  two  difficulties  were  encountered  in  the  operation  of 
the  tunnel.    When  velocities  above  6  miles  per  hour  were  introduced,  some 
deteriorating  plywood  sections  of  the  tunnel  began  to  disintegrate.  At 
lower  rates,  naturally  occurring  winds  introduced  significant  deviations 
to  the  wind  flows  in  the  tunnel.    This  latter  effect  could  be  minimized  by 
operating  at  night  when  the  air  was  relatively  calm. 

The  first  inspection  of  the  wind  tunnel  for  use  in  these  experiments 
was  undertaken  in  the  summer  of  1957.    A  preliminary  visit  was  planned  to 
check  on  the  conditions  of  the  equipment  and  on  the  feasibility  of  using 
the  fallout  disperser  and  burners  in  the  tunnel.    Because  of  the  excellent 
condition  in  which  the  tunnel  was  found,  sufficient  time  was  available  to 
proceed  with  a  number  of  extemporaneous  experiments.    Although  more  care- 
ful planning  would  have  simplified  subsequent  analysis  of  the  data,  it 
turned  out  that  the  results  of  these  preliminary  experiments  were  suffi- 
cient to  essentially  exhaust  the  conclusions  that  could  be  reached  with 
the  apparatus  available.    Thus,  no  further  trips  to  the  wind  tunnel  were 
necessary. 

Experimental  Procedures 

For  the  experimental  determination  of  the  influence  of  fire-induced 
convection  column  on  fallout  patterns,  the  flock  gun  was  mounted  at  the 
entrance  to  the  test  section  of  the  wind  tunnel,  centered  horizontally,  and 
6  inches  below  the  tunnel  ceiling.    The  burners  were  positioned  at  varying 
distances  ranging  from  3-1/2  feet  to  11  feet  downwind  of  the  flock  gun.  The 
floor  of  the  tunnel  was  covered  by  a  sheet  of  glossy  black  paper  30  inches 
wide  (preliminary  investigation  indicated  that  the  fallout  showed  up  much 
better  against  the  black  background  than  against  a  white  background). 
Filter  paper  collectors  6  inches  in  diameter  were  then  mounted  by  thumbtack 
along  the  center  of  the  black  paper  (and  center  of  the  wind  tunnel)  at  1- 
or  2-foot  intervals.    Additional  collectors  were  mounted  centered  10  inches 
to  the  right  and  left  of  the  central  row. 

The  experiment  then  consisted  of  adjusting  the  wind  tunnel  velocity 
to  2,  4,  or  6  miles  per  hour,  firing  the  flock  gun  containing  one  of  the 
five  size  fractions  of  fallout,  and  determining  the  resultant  pattern  by 
weighing  the  material  collected  on  the  filter  paper  collectors  to  the  near- 
est milligram  (samples  collected  ranged  from  2  or  3  grams  down  to  less  than 
a  milligram).    The  burners,  set  in  one  of  the  two  configurations  and  at  one 
of  four  distances  down  the  tunnel,  were  then  ignited  and  the  experiment 
repeated  with  a  convection  column  formed.    Observation  through  several  of 
the  optical  ports  in  the  walls  of  the  tunnel  confirmed  the  fact  that  the 
fallout  cloud  descended  until  it  reached  the  convection  column,  then  rose 
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sharply  and  finally  descended,  dispersed  over  a  larger  area  further  down- 
wind.   Typical  patterns  observed  following  runs  without  and  with  fire  are 
shown  in  figure  12.    A  limited  number  of  experiments  were  performed  with 
a  single  Meeker  burner  producing  a  much  smaller  flame.    The  resultant  change 
in  pattern  covered  too  small  an  area  for  detection  with  the  gross  measure- 
ments made  in  these  experiments.    However,  the  effect  of  this  small  burner 
is  illustrated  in  the  pattern  shown  in  figure  13. 

DISCUSSION 

SELECTION  OF  CHARACTERIZING  PARAMETERS 

Fallout  distribution  patterns  are  three-dimensional.    For  the  wind 
tunnel  experiments,  a  downwind  and  crosswind  dimension  serves  to  identify 
a  point  within  the  pattern,  and  the  amount  of  material  landing  at  this 
point  is  given  by  the  weight  of  sample  collected.    For  the  comparison  of 
patterns  produced  in  a  large  number  of  experiments,  a  point-by-point  anal- 
ysis, however,  is  not  feasible.     Instead,  each  pattern  must  be  characterized 
in  terms  of  a  limited  number  of  pertinent  parameters  and  these  parameters 
then  compared  as  the  experimental  conditions  are  changed.    For  purposes  of 
this  study,  for  example,  the  characterizing  parameters  can  be  compared  in 
the  presence  and  absence  of  fire,  all  other  factors  being  held  constant. 

Although  the  convection  columns  introduced  both  a  longitudinal 
(downwind)  and  transverse  (crosswind)  dispersion,  it  was  felt  that  the  pat- 
terns could  be  sufficiently  characterized  by  a  two-dimensional  measurement. 
Thus,  downwind  displacement  could  be  identified  by  downwind  measurements, 
and  additional  crosswind  dispersion  could  be  recognized  by  a  decrease  in 
the  total  amount  in  a  downwind  slice  through  the  pattern.    It  was  decided, 
therefore,  to  attempt  to  characterize  the  fallout  patterns  in  terms  of  the 
measurements  made  on  the  central  slice  taken  longitudinally  down  the  tunnel. 
Attempts  to  fit  this  downwind  distribution  in  terms  of  polynomial  expres- 
sions through  cubic  proved  unsuccessful.    The  data  fit  a  quadratic  expres- 
sion nicely  in  the  logarithm  of  the  amount  deposited.    However,  a  least 
squares  method  of  curve  fitting  (using  the  logarithms)  was  most  unsatis- 
factory; such  a  procedure  weighted  the  errors  in  measurements  of  the  very 
small  samples  at  the  tails  of  the  curve  equally  with  the  errors  in  the 
much  more  accurate  measurements  made  on  the  large  samples  near  the  peak. 

The  excellent  fit  to  a  quadratic  in  logarithms  led  to  the  estab- 
lishment of  the  characteristic  log  normal  distribution  to  identify  a  parti- 
cular pattern.    That  is,  it  was  found,that  the  weights  taken  in  the  central 
slice  fit  very  nicely  the  expression-' 

=  .4343(Ey)(Ax)  e  -  (log  x  -  log  u)2  (4) 
V2T  (a)(x)  2a^ 


2/    This  expression  can  be  transformed  simply  into  the  form 
log  y  =  a  +  b  (log  x)  +  c  (log  x)  . 
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Figure  12. --Fallout  patterns  with  and  without  fire.    Particle  size, 
105-149  H,  wind  velocity  2  m.p.h.    Left  Run  No.  25,  without 
fire;  right,  Run  No.  26,  with  fire. 
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in  which    y    is  the  weight  of  the  material  at  any  point  in  this  downwind 
slice;  Zy  ,  the  total  weight  of  the  samples  in  the  slice;    x  ,  the  down- 
wind distance;  and    Ax  ,  the  class  interval  (the  distance  between  points 
of  measurement).    The  characterizing  parameters,  then,  could  be    |i  ,  the 
logarithmic  mean  distance,  and    a  ,  the  standard  deviation  of  this  dis- 
tance, and  Zy  ,  the  total  amount  of  material  in  the  slice .  Alternative 
parameters  whose  concepts  are  more  readily  grasped  would  be    xp  ,  the 
distance  to  the  peak  in  the  parameter  curve,  which  is  related  to    |i  and 
a    by  the  expression 

log  x    =  log  u  -  2.3026a2  (5) 

and  yp  ,  the  amount  of  material  at  the  peak,  which  can  be  simply  derived 
from  equations  4  and  5  as 

v    -  .1733 far)  (Ax)  _  -  2.651a2 
yp  "   (a)(xp)  -e 

Thus,  the  downwind  displacement  of  the  peak  fallout  concentration  by  fire 
can  be  seen  directly  by  comparison  of  the  values  of    Xp  „    The  reduction 
in  amount  at  the  peak  can  be  seen  by  comparison  of  the  values  of    yp  . 
Of  this  reduction,  the  amount  resulting  from  downwind  dispersion  can  be 
seen  from  the  value  of    a  ,  the  remainder  being  due  to  crosswind  dis- 
persion. 

SAMPLE  ANALYSIS 

The  analytical  procedure  used  in  establishing  the  characteristic 
parameters  for  an  individual  run  has  considerable  applicability  in  many 
other  problems  and  will  be  described  in  detail  in  a  subsequent  report. 
As  applied  to  these  experiments  the  procedure  can  be  illustrated  briefly 
considering  one  of  the  more  than  a  hundred  experiments  performed.  The 
first  two  columns  of  table  2  present  the  basic  measurements  made  for 
experiment  number  108.    Column  1  gives  the  sample  distance  in  feet,  that 
is,  the  distance  between  the  fallout  disperser  and  the  midpoint  of  the 
fallout  collector  mounted  in  the  central  slice  through  the  tunnel. 
Column  2  gives  the  weight  of  fallout  collected  in  grams.    For  most  runs 
samples  were  taken  at  2-foot  intervals.     In  this  experiment,  to  more 
clearly  establish  the  position  of  the  peak,  additional  samples  were 
taken  at  1-foot  intervals  for  part  of  the  distance.    The  weights  given 
in  parenthesis  in  Column  2  represent  values  at  points  at  which  measure- 
ments were  not  made. 

To  establish  the  characterizing  parameters  for  this  run,  a  curve 
of    y    vs .     x    could  be  fitted  to  the  data  given  in  Columns  1  and  2  of 
table  2.    Since  such  a  curve  is  far  from  linear,  graphical  curve  fitting 
would  be  subject  to  considerable  error,  particularly  in  establishing  the 
position  of  the  peak.    However,  the  integral  of  a  log  normal  curve  may 
be  plotted  as  a  straight  line  on  log  probability  paper .    By  making 
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measurements  at  the  center  of  a  number  of  class  intervals  of  equal  size 
and  plotting  the  cumulative  values,,  expressed  in  percent,  against  the 
upper  class  limit  of  the  interval,  a  straight  line  should  result  on  log 
probability  paper.    Since  measurements  were  not  taken  at  equal  class  in- 
tervals, two  alternatives  are  available.    First,  only  measurements  made 
at  the  2-foot  intervals  could  be  used  (the  intervening  1-foot  measure- 
ments being  ignored),  and  the  result  obtained  in  this  way  is  given  in 
the  next  three  columns,  in  which  Column  3  gives  the  distance  to  the  upper 
class  limit,  Column  4  the  cumulative  weight  of  these  samples,  and  Column  5 
this  cumulative  weight  expressed  as  the  percentage  of  the  total  weight  of 
all  of  the  samples.    These  values  are  plotted  as  the  triangular  points  in 
figure  14. 

An  alternative  procedure  is  to  interpolate  values  for  the  1-foot 
intervals  at  which  measurements  were  not  made  and  plot  the  data  in  terms 
of  1-foot  class  limits.    The  values  given  in  parentheses  in  Column  2 
represent  linear  interpolations  between  the  adjoining  points  at  which 
measurements  were  made.    Column  6,  then,  gives  the  upper  class  limit  as- 
suming 1-foot  intervals  and  Columns  7  and  8  the  cumulative  weight  and 
cumulative  percent  for  this  case.    These  values  are  plotted  as  the  square 
points  in  figure  15. 

The  points  give  an  excellent  straight  line  in  the  central  portion 
of  the  graph,  but  turn  off  sharply  at  the  limits  (particularly  the  upper 
limit).    This  deviation  from  linearity  is  caused  by  the  fact  that  a  cer- 
tain amount  of  material  extends  beyond  the  limits  at  which  measurements 
were  made  (the  sum  of  all  the  measured  material  being  less  than  the  total 
area  under  such  a  distribution  curve) .    A  correction  may  be  made  for  this 
truncation  by  adding  amounts  for  distances  up  to  that  at  which  the  first 
measurement  was  made  and  beyond  that  at  which  the  last  measurement  was 
made.    This  has  been  done  in  Columns  9  and  10  and,  as  may  be  seen,  a 
straight  line  results  throughout  the  entire  range  of  measurements. 

Although  this  correction  for  truncation  is  arbitrary,  the  curve 
plotted  on  log  probability  paper  is  an  excellent  indication  of  the  cor- 
rectness of  the  values  chosen.     In  the  example  given,  a  correction  of 
about  1  percent  in  the  total  weight  was  necessary  to  straighten  out  the 
ends  of  the  curve.    An  overcorrection  of  the  same  amount  would  produce 
an  equal  curvature  in  the  opposite  direction.    Thus,  an  iterative  pro- 
cedure can  be  used  in  making  this  correction,  with  the  curve  plotted 
after  each  iteration  establishing  whether  an  additional  correction  is 
required.    The  accuracy  of  establishing  the  correct  curve  is  a  function 
of  the  degree  of  truncation;  however,  it  has  been  found  that,  as  long 
as  data  are  obtained  on  both  sides  of  the  peak,  a  remarkably  narrow 
range  of  alternative  curves  will  be  permitted  by  the  data.     Even  in  the 
highly  truncated  cases  where  measurements  are  made  on  only  one  side  of 
the  peak  (for  example,  in  these  experiments  when  most  of  the  fallout 
extended  beyond  the  limit  of  the  test  section  of  the  wind  tunnel)  fairly 
close  limits  can  be  set  upon  alternative  possibilities  by  considering 
such  factors  as  the  maximum  range  possible  in  values  of  the  cumulative 
weight . 
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Once  the  straight  line  given  in  figure  14  has  "been  established,  the 
characterizing  parameters  may  be  obtained  as  follows:    The  value  of    |i  may 
be  read  directly  as  the  value  of  the  ordinate  where  the  curve  crosses  the 
50-percentile  line.    The  value  of    a    is  obtained  from  the  slope  of  the 
curve.    For  example,  since  99  percent  of  the  material  is  found  within  the 
limits    ±  2.576a    of  the  logarithmic  mean,     a    can  be  obtained  by  taking 
the  difference  between  the  logarithm  of  the  ordinate  as  the  curve  crosses 
the  99.5  percent  line  and  as  the  curve  crosses  the  0.5  percent  line  and 
dividing  by  5.152.    The  position  of  the  peak,  then,  can  be  obtained  from 
the  values  of    \x    and    a    using  equation  5.    It  is  interesting  to  note 
that  in  establishing  the  best  straight  line  through  the  points,  a  line 
leading  to  a  higher  value  of    \i    invariably  has  also  a  steeper  slope. 
Thus,  the  value  of    Xp    may  be  quite  firmly  established  with  larger  errors 
in  both    |j    and    a  .    Having  established  the  value  of    a    and    |_i  ,  a  curve 
of    y  vs.  x    may  be  plotted  directly  from  equation  4»    Such  a  curve  is 
given  in  figure  15,  together  with  the  original  data  shown  in  Columns  1  and 
2  of  table  2.     It  may  be  seen  that  the  technique  of  obtaining  the  charac- 
teristic parameters  from  the  integrated  curve  can  smooth  out  rather  size- 
able point-to-point  variations  in  the  differentiated  curve  and  permit  the 
accurate  establishment  of  the  position  of  the  peak  when  a  freehand  graph- 
ical fitting  would  be  subject  to  considerable  error.    The  value  of  yp 
may,  of  course,  be  obtained  from  figure  15  or  may  be  computed  directly 
from  equation  6. 

EXPERIMENTAL  RESULTS 

The  values  of  the  characterizing  parameters  of  all  the  experiments 
are  presented  in  table  3.     Inspection  of  these  values  shows  the  obvious 
effect  of  fire,  of  particle  size,  and  of  wind  velocity.    The  lighter, 
slower -falling  particles  travelled  farther  down  the  tunnel  before  land- 
ing.   The  stronger  the  wind,  the  farther  down  the  tunnel  the  particles 
went.    In  practically  every  case,  the  effect  of  the  fire  was  to  introduce 
a  further  lateral  dispersion  and  downwind  movement  of  the  fallout,  giving 
a  more  disperse,  less  concentrated  pattern.    The  one  obvious  exception  to 
this  effect  of  fire  occurred  in  the  few  experiments  in  which  the  burner 
was  located  at  the  edge  of  the  no- fire  pattern.     In  these  cases  air  flow 
into  the  convection  column  shifted  the  peak  of  the  pattern  in  an  upwind 
direction. 

Cursory  inspection  of  the  numbers  in  table  3  did  not  show  whether 
there  was  any  significant  effect  of  burner  position  at  the  3-1/2  and  7- 
foot  positions  or  of  burner  configuration  (in  the  parallel  vs.  V-shaped 
configurations).    Consequently,  a  statistical  analysis  was  undertaken  to 
determine  whether  these  factors  introduced  a  significant  effect.  The 
experiment  was  considered  to  be  a  five-factor  factorial.    Eliminating  all 
levels  with  insufficient  replicates,  the  analysis  was  confined  to  two 
levels  of  each  of  four  factors  (burner  position,  burner  configuration, 
wind,  and  fire  condition)  and  four  levels  of  the  remaining  factor  (par- 
ticle size).    Owing  to  the  extreme  disproportionality  of  frequencies  of 
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replicates,  it  was  decided  that  a  complete  least-square  solution  should  be 
attempted.    Out  of  the  possible  63  constants  in  an  experiment  such  as  this, 
the  large  number  of  cells  without  observations  reduced  the  possible  number 
of  constants  to  44-    Again,  owing  to  the  way  that  observations  were  distri- 
buted among  the  possible  factor  combinations,  models  containing  certain 
constants  resulted  in  indeterminant  solutions  (singular  matrices).  Thus, 
no  models  were  possible  containing  interactions  involving  burner  condition, 
nor  containing  interactions  greater  than  first  order  with  burner  configu- 
ration. 

For  the  seven  models  that  were  admissible,  multiple  correlation  co- 
efficients were  computed  for    Xp,  \±,  and  a,  using  the  IBM  701  machine  avail- 
able at  the  University  of  California  Computer  Center.     On  the  basis  of  these 
results  (table  4) >  the  model  containing  particle  size,  wind  velocity,  and 
fire  condition  alone  (including  interactions)  was  adopted.    Thus,  in  obtain- 
ing average  values  of  the  characteristic  parameters,  the  burner  configura- 
tion and  burner  position  (for  those  runs  in  which  the  burner  was  at  3-1/2 
or  7  feet)  were  ignored. 

The  average  values  of  the  characteristic  parameters  were  then  obtained 
for  a  given  fire,  wind,  and  particle-size  situation  (table  5).  Examination 
of  the  values  of    Xp  and  yp    in  this  summary  table  show  clearly  the  effect  of 
the  fire-induced  convection  column  in  these  experiments.     In  each  case,  fire 
resulted  in  a  marked  downwind  displacement  of  the  peak  in  the  fallout  distri- 
bution pattern  (and  thus  of  the  whole  pattern)  with  a  decrease  in  value  at 
the  peak  ranging  up  to  a  factor  of  10.     How  this  effect  was  influenced  by 
the  walls  and  ceiling  of  the  wind  tunnel  is  not  clear.    Consequently,  these 
results  are  not  directly  extrapolatible  to  free-air  conditions. 

One  final  measurement  bearing  upon  the  interpretation  of  these  ex- 
perimental runs  is  worthy  of  mention.    Part  of  the  downwind  distribution  of 
fallout  in  any  run  may  be  attributed  to  further  fractionation  of  the  parti- 
cles within  the  particular  sieve  size  fraction  on  the  basis  of  fall  velocity. 
Consequently,  for  a  number  of  runs  samples  were  collected  on  transparent 
gummed  paper  mounted  along  the  centers  of  the  tunnel.    Representative  samples 
were  then  measured  under  the  microscope  using  the  statistical  technique  de- 
scribed previously.     The  results,  as  would  be  expected,  show  a  definite 
fractionation  down  the  tunnel,  with  mean  particle  size  decreasing  with  in- 
creasing distance  downwind  (table  5). 
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Table  4 • --Multiple  correlation  coefficients 


for  admissible  models 


Number 
of 

: Multiple 

correlation 
for 

coefficient 

Model  : 

variables 

:  |i 

:  a 

1st  order ,  including: 

rire,  winQ,  parxicie 
size 

D 

VU.L) 

Aq  n 
Oy .  U 

77*  3 

Fire,  wind,  size,  burner 
position 

7 

90.0 

89.2 

78.8 

Fire,  wind,  size,  burner 
configuration 

7 

90.4 

89.3 

77.3 

Fire,  wind,  size, 

position,  configuration 

8 

90.4 

89.5 

78.8 

2nd  order  (1st  order 
interactions ) : 

Fire,  wind,  particle 
size 

13 

92.5 

91.7 

80.1 

Fire,  wind,  size, 
configuration 

19 

93.8 

92.8 

81.5 

3rd  order  (2nd  order 
interactions) : 

Fire,  wind,  size 

16 

94.6 

94.3 

82.3 
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CONCLUSION 


It  has  been  shown  that,  under  some  conditions  with  sufficiently  high 
updraft,  a  fire-induced  convection  column  can  have  a  marked  effect  upon 
radiological  fallout  patterns.    The  effect  of  the  column  itself  will  "be 
markedly  influenced  by  the  air  flow  in  and  around  the  column.    The  experi- 
ments conducted  in  this  study  were  performed  in  a  wind  tunnel  in  which  this 
flow  was  influenced  by  the  tunnel  walls  and  ceiling.    Consequently  the  re- 
sults are  not  directly  applicable  to  open-air  conditions.    However,  the 
results  do  suggest  quite  strongly  that  a  marked  effect  will  be  found  for 
a  large-scale  fire  burning  freely.    Under  conditions  of  sizeable  civil  de- 
fense activity,  such  an  effect  could  greatly  alter  the  necessary  civil 
defense  planning  and  action.    Consequently,  open-air  experiments  on  a 
larger  scale  should  be  performed  in  an  attempt  to  establish  the  parameters 
affecting  the  fallout  pattern  in  the  presence  of  a  fire  in  the  event  of 
a  nuclear  detonation  where  fallout  occurs  and  combustible  material  is 
available. 
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